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ABSTRACT 

Aims. We present basic atmospheric parameters (T eS , log g, y t , and [Fe/H]) as well as luminosities, masses, radii, and absolute radial 
velocities for 348 stars, presumably giants, from the ~1000 star sample observed within the Penn State-Torun Centre for Astronomy 
Planet Search with the High Resolution Spectrograph of the 9.2 m Hobby-Eberly Telescope. The stellar parameters (luminosities, 
masses, radii) are key to properly interpreting newly discovered low-mass companions, while a systematic study of the complete 
sample will create a basis for future statistical considerations concerning the appearance of low-mass companions around evolved 
low- and intermediate-mass stars. 

Methods. The atmospheric parameters were derived using a strictly spectroscopic method based on the LTE analysis of equivalent 
widths of Fe I and Fe II lines. With existing photometric data and the HlPPARCOS parallaxes, we estimated stellar masses and ages 
via evolutionary tracks fitting. The stellar radii were calculated from either estimated masses and the spectroscopic log g or from 
the spectroscopic T e ff and estimated luminosities. The absolute radial velocities were obtained by cross-correlating spectra with a 
numerical template. 

Results. We completed the spectroscopic analysis for 332 stars, 327 of which were found to be giants. A simplified analysis was 
applied to the remaining 16 stars, which had incomplete data. The results show that our sample is composed of stars with effective 
temperatures ranging from 4055 K to 6239 K, with log g between 1.39 and 4.78 (5 dwarfs were identified). The estimated luminosities 
are between log L/L e = -1.0 and 3 and lead to masses ranging from 0.6 to 3.4 M e . Only 63 stars with masses larger than 2 M e were 
found. The radii of our stars range from 0.6 to 52 R a with the vast majority between 9-11 R a . The stars in our sample are generally 
less metal-abundant than the Sun with median [Fe/H] = -0.15. The estimated uncertainties in the atmospheric parameters were found 
to be comparable to those reached in other studies. However, due to lack of precise parallaxes, the stellar luminosities and, in turn, 
the masses are far less precise, within 0.2 M e in best cases and 0.3 M e on average. 

Key words, stars: fundamental parameters - stars: atmospheres - stars: late-type - stars: abundances - planetary systems 


1. Introduction 

Since the discovery of the first extrasolar planetary systems by 
Wolszczan & Frail (1992), Mayor & Queloz (1995), and Marcy 
& Butler (1996), over 700 planets have been found around other 
stars. The richness of exoplanets and the architectures of plan¬ 
etary systems emerging today are amazing and raise questions 
about a general picture of planet formation and evolution. But 
before such a general picture becomes available, studies of plan¬ 
etary systems in various environments are essential. 

The observational techniques applied in searches for exo¬ 
planets are all limited to their areas of competence, i.e., the 
range in stellar or planetary parameters (effective temperatures, 


* Based on observations obtained with the Hobby-Eberly Telescope, 
which is a joint project of the University of Texas at Austin, 
Pennsylvania State University, Stanford University, Ludwig- 
Maximilians-Universitat Miinchen, and Georg-August-Universitat 
Gottingen. 

** Tables 1 and 5 are only available at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) orvia 
http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/547/A91 


masses, planetary periods or orbital separations, etc.) that are 
available to them. Due to the limitations of current detectors, 
the direct imaging searches for exoplanets typically deliver 
massive planets at large orbital distances from low-mass stars 
(cf. 2M1207 b - Chauvin et al. 2004; or HD 8799 system - 
Marois et al. 2008). The transit planet searches, conveniently 
operating in the sin i ~ 1 domain (i being orbital inclination), 
discover exoplanets with relatively short orbital periods (orbital 
separations) due to the natural limit of the observable transit 
length set by the lengths of a night (a condition now being re¬ 
laxed with the multi-site or space observations, cf. HD 80606 b - 
Hidas et al. 2010; Hebrard et al. 2010; Kepler 1 lg — Lissauer 
et al. 2011). Searches for planets based on the radial velocity 
(RV) technique are the most versatile. They detect planetary can¬ 
didates at various orbital separations and in a wide range of plan¬ 
etary minimum masses, as long as sin/ + 0 (cf. GJ 581 sys¬ 
tem - Mayor et al. 2009; and 47 UMa d - Gregory & Fischer 
2010). Even this method has its limitations, though. Since it 
requires a large number of narrow spectral lines, it focuses on 
the slowly rotating cool stars of late-F to M spectral type which 
sets an upper limit to the main sequence planetary hosts masses 
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(M < 1.5 M 0 ). Unfortunately this technique is also seriously af¬ 
fected by the stellar activity (Queloz et al. 2001), which, if not re¬ 
solved, adds noise (jitter) and limits the precision in RV (Carney 
et al. 2003; Hekker et al. 2006; Cumming et al. 2008). It is, how¬ 
ever, the most suitable in searches for planets around stars more 
massive than ~1.5 M Q . In terms of evolution, these stars slow 
down their rotation and decrease their effective temperatures, be¬ 
coming red giants available for precise RV measurements. 

The Penn State-Torun Centre for Astronomy Planet Search 
(PTPS) is aiming to detect and characterize planetary sys¬ 
tems around stars at various evolutionary stages (Niedzielski 
& Wolszczan 2008) with the long-term goal of describing the 
evolution of planetary systems with aging stars. To this end, 
over 1000 stars are monitored with the Hobby-Eberly Telescope 
(HET) for RV variations using the high-precision iodine-cell 
technique (Butler et al. 1996). The sample is mainly composed 
of evolved low-mass and intermediate-mass stars: —575 giants 
(including 348 stars from the red giant clump (RGC) sample 
presented here) and -225 sub-giants, but it also contains -200 
slightly evolved dwarfs. Several stars with planetary-mass com¬ 
panions have already been discovered, mainly around red gi¬ 
ants (Niedzielski et al. 2007, 2009a,b; Gettel et al. 2012b, a). 
These findings supplement the slowly growing population of 
-50 planets currently known around evolved stars, discov¬ 
ered in such projects as the McDonald Observatory Planet 
Search (Cochran & Hatzes 1993; Hatzes & Cochran 1993), the 
Okayama Planet Search (Sato et al. 2003), the Tautenburg Planet 
Search (Hatzes et al. 2005), the Lick K-giant Survey (Frink et al. 
2002), the ESO FEROS planet search (Setiawan et al. 2003a, b), 
the Retired A Stars and Their Companions (Johnson et al. 2007), 
the CORALIE and HARPS search (Lovis & Mayor 2007), and 
the Boyunsen Planet Search (Lee et al. 201 1). 

With this paper, we are starting a series devoted to a detailed 
description of stars incorporated in the complete PTPS sample, 
their atmospheric parameters, elemental abundances, rotation 
velocities, kinematical properties, masses, radii, ages, etc. The 
derived parameters not only describe our sample but also allow 
for proper interpretation of the planet search results. The stellar 
masses are essential in extracting minimum masses of the plan¬ 
etary mass companions (wp sin /), while the radii and projected 
rotation velocities are key ingredients in stellar activity consider¬ 
ations etc. Since the complete sample contains -1000 stars ran¬ 
domly distributed in the skies, it is of interest for general stellar 
and galactic astrophysics as well. It may be used for studies of 
galactic structure and evolution, galactic distribution of plane¬ 
tary systems (Haywood 2009), and the galactic habitable zone 
(Gonzalez et al. 2001 ; Lineweaver et al. 2004; McCabe & Lucas 
2010; Gowanlock et al. 201 1). 

The paper aims to determine atmospheric parameters, such 
as effective temperatures (T e g), stellar gravitational acceler¬ 
ations (log g), microturbulence velocities (p t ), and metallici- 
ties ([Fe/H]) for 348 GK-type stars, presumably red clump giants 
observed within the PTPS survey. Together with existing photo¬ 
metric data and parallaxes, they will allow us to estimate stellar 
masses (M/M Q ), radii (R/R e ), and ages. 

The RGC, together with the main sequence (MS) and the red 
giant branch (RGB), is one of the most characteristic features 
of the Hertzsprung-Russell diagram (HRD). It was identified by 
Cannon (1970) and Faulkner & Cannon (1973) as the location of 
steady-core helium burning. That finding was elaborated in de¬ 
tail by Girardi (1999). Jimenez et al. (1998) used the RGC stars 
to estimate the age of the Galaxy. Due to very similar intrinsic 
brightness of all stars at that stage, the RGC was proposed as a 
standard candle in distance estimates and applied to the Galactic 


center (Paczynski & Stanek 1998), the Large Magellanic Cloud 
(Salaris et al. 2003), and other stellar systems. 

Their interesting stellar evolution phase and general as- 
trophysical application make the red giants, and especially 
the RGC stars, the subject of numerous surveys. McWilliam 
(1990) obtained T e s for 671 stars from broadband photometry, 
log g via isochrone fitting, and the local thermodynamical equi¬ 
librium (LTE) abundances from high-resolution spectra. Zhao 
et al. (2001) delivered atmospheric parameters, iron abundances, 
a-element enhancements, and masses for 39 red clump giants 
selected from the HlPPARCOS catalogue. Famaey et al. (2005) 
provided the CORAVEL radial velocities for a sample of about 
6600 K-type giants. Detailed spectroscopic studies of various 
subsamples of the galactic red clump giants were presented in 
Mishenina et al. (2006), Bizyaev et al. (2006, 2010), Hekker 
& Melendez (2007), Liu et al. (2007), Luck & Heiter (2007), 
Takeda et al. (2008), and Puzeras et al. (2010). 

Gontcharov (2008) summarized kinematical properties and 
galactic distribution of 97 348 RGC stars based on Tycho-2 
proper motions and Tycho-2 and 2MASS photometry. Valentini 
& Munari (2010) studied 277 RGC stars and presented accu¬ 
rate multi-epoch radial velocities, atmospheric parameters, dis¬ 
tances, and space velocities ({/, V, W ). Moreover, Tautvaisiene 
et al. (2010) in a recent paper presented 12 C/ 13 C abundance ratio 
determinations in a sample of 34 Galactic clump giants as well 
as abundances of nitrogen, carbon, and oxygen. 

Meanwhile, new insight into giants, including clump giants, 
was gained through CoRoT (Baglin et al. 2006) and Kepler 
(Gilliland et al. 2010) precise photometry. The solar-like oscil¬ 
lations in red giants were first identified in ground-based RV ob¬ 
servations (Hatzes & Cochran 1994; Merline 1999). De Ridder 
et al. (2009) reported the presence of radial and non-radial oscil¬ 
lations in over 300 giants observed with CoRoT. These oscilla¬ 
tions were used in massive masses and radii determinations for 
giants (Kallinger et al. 2010; Bedding et al. 2010; Hekker et al. 
2011 ). 

In Sect. 2 of this paper, we describe the sample and the ob¬ 
servational material. The procedure of Fe lines equivalent widths 
measurements and its precision are presented in Sect. 3. We give 
the outline of the method applied in our spectroscopic analysis 
in Sect. 4, while in Sect. 5, we present several tests to evaluate 
the reliability of the method and estimate uncertainties. Section 6 
contains a summary of the atmospheric parameters determined 
for our stars. These parameters are used in Sect. 7 to estimate 
stellar integrated parameters: luminosities, masses, radii, and 
ages. Section 8 presents the method and results of our RV de¬ 
terminations, while a comparison of our results with other pub¬ 
lished data, is found in Sect. 9. Section 10 discusses atmospheric 
parameters, and conclusions are presented in Sect. 11. 


2. Target selection and observations 

The sample of 348 stars presented here constitute the RGC sub¬ 
sample of the PTPS. Originally the SIM-EPIcS 1 reference can¬ 
didates, they were selected as presumable RGC stars due to their 
photometry and reduced proper motions (Gelino et al. 2005 ; Law 
et al. 2006). 


1 Space Interferometry Mission (SIM) key project: Extra-solar Planet 

Interferometric Survey (EPIcS) 

http://planetquest.jpl.nasa.gov/SIM 
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Fig. 1. Histograms of the apparent magnitudes in V band (left panel) 
and (B - V) (right panel) for the 348 red clump giants observed within 
the PTPS. 


They are relatively bright held stars, with V between 4.5 and 
11 mag, but in the SIMBAD 2 database, spectral types and/or lu¬ 
minosity class estimates are available only for a few (four are 
classified as dwarfs and 22 as giants). The basic photometric data 
for the sample stars come from the 2MASS Point Source cata¬ 
log (Cutri et al. 2003) and the Tycho-2 catalog (H0g & Murdin 
2000), compiled by Gelino et al. (2005). For several stars, miss¬ 
ing data were gathered from Kharchenko & Roeser (2009). The 
distributions of visual magnitudes and B - V color indices are 
presented in Fig. 1. The vast majority, 92% of stars, are fainter 
than V = 8 mag, with (B— V) between 0.7 and 1.8 mag (1.1 mag 
on average). 

Our sample is composed of much fainter stars than those in¬ 
cluded in most of the other RGC surveys (Famaey et al. 2005; 
Mishenina et al. 2006; Hekker & Melendez 2007; Takeda et al. 
2008; Valentini & Munari 2010) and can be compared to the 
more numerous sample of Bizyaev et al. (2006, 2010). The list 
of stars, their identification, and basic observational parameters 
are presented in Table 1. 

Our high-quality, high-resolution optical spectra have been 
collected since 2004 with the HET (Ramsey et al. 1998), located 
in the McDonald Observatory. The telescope was equipped with 
the High Resolution Spectrograph (HRS) fed with a 2 arcsec 
fiber, working in the R — 60000 resolution (Tull 1998). The 
observations were performed in the queue-scheduled mode 
(Shetrone et al. 2007). The spectra consisted of 46 echelle or¬ 
ders recorded on the “blue” CCD chip (407-592 nm) and 24 or¬ 
ders on the “red” one (602-784 nm). The signal-to-noise ratio 
was typically better than 200-250 per resolution element. Basic 
data reduction (flat fielding, extraction of the spectrum, wave¬ 
length calibration, normalization to continuum) was performed 
using standard IRAF 3 tasks and scripts. The wavelength calibra¬ 
tion was based on Th-Ar lamp spectrum identification presented 
in Mader & Shetrone (2002). 

3. Equivalent widths of Fe lines 

Since equivalent widths (EWs) of hundreds of iron lines 
for every object were needed, we used the DAOSPEC 4 

2 The SIMBAD astronomical database is operated at CDS, Strasbourg, 
France. 

3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

4 http://www2.cade-ccda.hia-iha.nrc-enrc.gc.ca/ 
community/STETSON/daospec 

http://www.bo.astro.it/~pancino/proj ects/daospec.html 



EW DA ospec [mA] 

Fig. 2. Comparison of the automatic (DAOSPEC) and handmade 
(IRAF) EWs measurements of 5695 Fe lines for 29 stars. The solid 
line presents the one-to-one relation, the dashed one is the best linear 
fit to the data. The individual uncertainties of the DAOSPEC EWs mea¬ 
surements are presented. 

(Stetson & Pancino 2008) to measure the EWs automatically. 
However, to validate reliability of the results and estimate 
uncertainties in EWs, we measured the spectra of 29 stars 
manually with the IRAF (Gaussian fitting) as well. Figure 2 
compares the IRAF (manual) with the DAOSPEC (automatic) 
measurements. We found that the best agreement between the 
measurements is reached for lines with EWs lower than 200 mA. 
Lines with larger EWs revealed some discrepancy towards 
higher manual EWs values. In general, the linear fit to all 
measurements is nearly a one-to-one relation. For 5632 lines 
with the EW < 200 mA, we found TTViraf = (0.937 ± 
0.005)TTV) JA osphc + (8.46 ± 0.35) and the Pearson correlation 
coefficient of R = 0.939. The mean difference between the 
IRAF and the DAOSPEC measurements was 5.3 mA with a scat¬ 
ter of 12.8 mA. The mean uncertainty in the DAOSPEC EWs 
was 2.13 mA. We proceeded with the automatic measurements 
in the further analysis. 

To apply the procedure described in the following section, 
we measured the EWs of neutral (Fel) and ionized (Fell) iron 
absorption lines for which precise atomic data, i.e., labora¬ 
tory wavelengths, lower excitation potentials (x ), and oscilla¬ 
tor strength (gf) were available from Takeda et al. (2005a,b). 
All these atomic data were compiled from the line-lists of 
Grevesse & Sauval (1999), Meylan et al. (1993), and Kurucz 
et al. (1984). Due to differences in the spectral coverage 
and CCD chip cosmetics, we were able to identify up to 
296 lines (269 Fel and 27 Fell lines) from the initial sam¬ 
ple of ~330 iron lines from Takeda et al. (2005a) in the 
HET/HRS spectra. After rejecting blends as well as very 
weak and saturated lines, we used the EWs of 190 Fel and 
15 Fe II lines per star on average. Only unblended lines with the 
EWs ranging from 5 mA to 200 mA were included in the final 
spectroscopic analysis. 

4. An outline of the analysis method 

The iron EWs were analyzed with the TGVIT (Takeda et al. 
2002a, 2005a), which is a part of the SPTOOL 5 package. This 
numerical algorithm is a modified version of an earlier code 

5 http://optik2.mtk.nao.ac.jp/~takeda/sptool 
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Fig. 3. Illustration of the TGVIT analysis of TYC 3018-00996-1: 
Fe abundance vs. x (top panel ) relation and Fe abundance vs. EW ( bot¬ 
tom panel ) relation. The open and filled circles stand for 176 Fel and 
12 Fell lines respectively. The mean Fe abundance is shown on each 
panel as a dashed line. 


if no relation between Fe abundance and these two parame¬ 
ters was present. Other conditions evaluated for every star were 
the equivalence of mean Fe abundances for Fe I and Fe II lines 
and the scatter around average abundance. After visual inspec¬ 
tion, the Fe lines showing apparent deviations from the mean 
Fe abundances were rejected. In the following step, the param¬ 
eters were redetermined and the procedure was repeated un¬ 
til coherent solutions for all parameters were reached. The re¬ 
jected lines were typically unrecognized blends or lines with 
EWs larger than 150 mA and the excitation potential of less 
than 0.5 eV. 

The Fe abundance depends on both T e t j and log g. The cou¬ 
pling is so strong that in the iterative approach of Takeda et al. 
(2002a) T e ff and log g are actually simultaneously fitted and 
treated as “one variable” ( v t being the other). This coupling may, 
in principle, lead to degeneracy of purely spectroscopic parame¬ 
ters because higher effective temperatures may be compensated 
by gravitational acceleration. However, in our sample, which is 
composed of stars with effective temperatures similar or lower 
to the Sun, we do not expect this effect to cause significant un¬ 
certainties. We refer the reader to Takeda et al. (2002a) for a 
detailed discussion of this problem and for an extensive discus¬ 
sion on the behavior of the D 2 function, along with its related 
components, around the converged solutions. 


and is based on the atmospheric models computed by Kurucz 
(1993a,b). The TGVIT uses an iterative procedure to achieve 
more cohesive parameter fits, and the version available to us was 
more suitable for lower effective temperatures (7’eii > 4250 K) 
than the original version (Takeda et al. 2002a). As a result, we 
could apply it to the late-type stars in our sample. More de¬ 
tails concerning the iteration procedures, data calculations, and 
formulation of the numerical problems are presented in Takeda 
et al. (2002a,b, 2005a, b). 

This purely spectroscopic method is based on an analysis of 
Fe I and Fe II lines and fulfills three conditions resulting from the 
assumption of the local thermodynamical equilibrium (LTE): 

1. the abundances derived from Fel lines may not show any 
dependence on the lower excitation potential x (excitation 
equilibrium); 

2. the averaged abundances from Fe I and Fe II lines must be 
equal (ionization equilibrium); 

3. the abundances derived from Fe I lines may not show any de¬ 
pendence on the EWs (matching the curve of growth shape). 

We used the TGVIT to calculate T e g, log g, v t , and [Fe/H], 
defined in the standard manner as [Fe/H] = \og(Nf e /Nn) - 
\og(Nf e /Nn)cs, with the solar value of 7.50 (Kurucz 1993a; 
Holweger et al. 1991). 

Because of the relatively low T e g of our stars, the number 
of available Fell lines was not sufficiently large (less than 10% 
of Fel lines). Consequently, we neglected the dispersion of 
Fe II abundances while estimating the minimization function D 2 
(TGVIT defines the dispersion D 2 as a function of three ar¬ 
guments: T e ff, log g, and u t that is minimized to obtain proper 
Fe abundance). For the coefficients included in D 2 , we assumed 
the values of c\ - 0, C 2 = 1, C 3 = 0 (see Takeda et al. 2002a, 
for details). 

In the vast majority of cases, all three above-mentioned con¬ 
ditions were fulfilled. This was checked for every star by visual 
evaluation of relations between Fe abundances and x or EWs 
like those presented in Fig. 3. The assumptions were justified 


5. Verification of the adopted approach 

Before running the TGVIT for our program stars, we performed 
a series of tests using various input data. First of all, we tested 
our installation of the TGVIT by running it with the data for the 
Sun and r Cet (HD 10700) distributed together with the code. 
These EWs were obtained from high-dispersion spectra obtained 
at the Okayama Astrophysical Observatory with the 188 cm tele¬ 
scope (Takeda et al. 2005a). Our results are compared with those 
of Takeda et al. (2005a) in Table 2. The results are practically 
identical, which proves that our installation works well. 

In the next step, we determined atmospheric parameters 
(i) for several stars, for which similar determinations were avail¬ 
able from Katz et al. (1998), using Echelle spectra available from 
the ELODIE Library (Soubiran et al. 1998) and (ii) for a few 
well-known stars with planets, for which recent atmospheric pa¬ 
rameter determinations are available in Butler et al. (2006), us¬ 
ing our own HET/HRS spectra. The purpose of these tests was 
to estimate reliability of the adopted method and to search for 
systematic effects. 

Another important reason for that analysis was to assess the 
reliability of uncertainties delivered by the TGVIT as these are 
determined from the quality of fits for T e ff, log g, and microtur¬ 
bulence velocity and should be considered intrinsic. Only the 
uncertainties in [Fe/H] are based on actual scatter between in¬ 
dividual determinations for all lines in use and are expected 
to be realistic. 

5. 1. Tests with ELODIE data 

With the TGVIT, we analyzed the spectra of 13 randomly se¬ 
lected stars, downloaded from the ELODIE Library (Soubiran 
et al. 1998) 6 at the Observatoire de Haute-Provence. Both giants 


6 Library of the ELODIE spectra of F5 to K7 stars is only available in 
electronic form via 

http://vizier.cfa.harvard.edu/viz-bin/VizieR? 
-source=J/A+AS/133/221 
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Table 2. Atmospheric parameters of the Sun and r Cet determined by Takeda et al. (2005a) and our results. 





Takeda et al. 



Our results 


Name 

Sp. 

T e g 

log g 

v t 

[Fe/H] 

Teff 

logg 

Vt 

[Fe/H] 


type 

[K] 

[cm s~ 2 ] 

[km s 1 ] 


[K] 

[cm s~ 2 ] 

[km s 1 ] 


Sun 

G2V 

5785 ± 8 

4.44 ± 0.02 

0.96 ± 0.06 

0.01 ± 0.04 

5785 ± 8 

4.44 ± 0.02 

0.96 ± 0.06 

0.02 ± 0.04 

r Cet (HD 10700) 

G8V 

5420 ± 0 

4.68 ± 0.00 

0.66 ±0.10 

-0.43 ± 0.06 

5420 ± 0 

4.68 ± 0.00 

0.65 ±0.10 

-0.43 ± 0.06 


Table 3. Atmospheric parameters determined by Soubiran et al. (1998) and our results for 13 stars observed with ELODIE. 


Name 

Sp. 

type 

Teff 

[K] 

Soubiran et al. 

log g 

[cm s~ 2 ] 

[Fe/H] 

T e ff 

[K] 


Our results 
I°g9 
[cm s~ 2 ] 

[Fe/H] 

HD 

10380 

K3III 

4057 ± 74 

1.43 ±0.25 

-0.25 ±0.14 

4192 

+ 

13 

2.03 

±0.06 

-0.36 ± 0.07 

HD 

10700 

G8V 

5264 ± 74 

4.36 ± 0.25 

-0.50 ±0.14 

5311 

± 

28 

4.40 

±0.07 

-0.41 ±0.10 

HD 

22049 

K2V 

5052 ± 74 

4.57 ± 0.25 

-0.15 ±0.14 

5113 

+ 

40 

4.65 

±0.10 

0.04 ±0.12 

HD 

30834 

K3III 

4115 ± 74 

1.73 ± 0.25 

-0.21 ±0.14 

4246 

+ 

15 

1.99 

±0.08 

-0.45 ± 0.08 

HD 

34411 

GOV 

5835 ± 63 

4.17 ± 0.20 

0.06 ±0.14 

5886 

± 

15 

4.20 

±0.03 

0.15 ± 0.05 

HD 

39003 

K0III 

4610 ± 74 

2.18 ± 0.25 

-0.10 ±0.14 

4618 

± 

30 

2.18 

± 0.11 

-0.03 ±0.14 

HD 

39853 

K3 

3994 ± 74 

1.00 ±0.25 

-0.40 ±0.14 

4094 

+ 

13 

1.63 

±0.06 

-1.00 ±0.08 

HD 

113226 

G8III 

4983 ± 74 

2.80 ± 0.25 

0.05 ±0.14 

5016 

+ 

28 

2.49 

±0.09 

0.11 ± 0.13 

HD 

124897 

K2III 

4361± 74 

1.93 ±0.25 

-0.53 ±0.14 

4305 

+ 

15 

1.77 

±0.07 

-0.55 ± 0.09 

HD 

132142 

K1V 

5108 ± 74 

4.50 ± 0.25 

-0.55 ±0.14 

5170 

± 

15 

4.64 

±0.04 

-0.35 ± 0.07 

HD 

175305 

G5III 

4899 ± 131 

2.30 ± 0.43 

-1.43 ±0.20 

5074 

+ 

30 

2.89 

±0.16 

-1.31 ±0.12 

HD 

204613 

G0III 

5727 ± 63 

3.84 ± 0.20 

-0.51 ±0.14 

5842 

+ 

18 

4.13 

±0.04 

-0.21 ±0.06 

HD 

217014 

G5V 

5729 ± 63 

4.12 ± 0.20 

0.11 ± 0.14 

5789 

± 

18 

4.34 

±0.05 

0.24 ± 0.07 



Fig. 4. Comparison of the atmospheric parameters for 13 stars observed with the ELODIE. The panels present the relations (top) and differ¬ 
ences ( bottom ) between our results and those obtained by Soubiran et al. (1998) for T t g (a), log g (b), and [Fe/H] (c). The solid lines represent 
one-to-one relations. 


and dwarfs with spectral types of G0-K3 were selected to cover 
parameter space similar to our sample. 

These spectra were obtained with a resolution of R - 42 000 
(Katz et al. 1998), which was lower than ours. The spectral 
wavelength range covered was 390-680 nm in 67 Echelle or¬ 
ders. Typical signal-to-noise ratio was 100-150 per pixel. The 
spectra were reduced, i.e., straightened, wavelength-calibrated, 
and cleaned from cosmic ray hits, bad pixels, and telluric lines. 
Because of the narrower spectral range, we were able to mea¬ 
sure 150 Fel and 15 Fell lines per star on average for the 
further analysis. 

The results obtained with the TGVIT for these 13 stars are 
presented and compared with results of Soubiran et al. (1998) in 
Table 3 and Fig. 4. 


The mean intrinsic uncertainties of our determinations are 
crT e ff = 21 K, crlog g - 0.07, and cr[Fe/H] = 0.09. According 
to Katz et al. (1998), the average uncertainties based on the 
internal accuracy assumed for stars with a signal-to-noise ra¬ 
tio around 100 are crT e g - 86 K, crlog g = 0.28, and 
cr[Fe/H] = 0.16. We note that the mean intrinsic uncertainties 
in T e g and log g as given by the TGVIT are about four times 
smaller than uncertainties estimated by Katz et al. (1998), while 
in the case of [Fe/H] they are comparable. An overall agreement 
of the results is clear from Fig. 4. 

We calculated the mean differences and scatter (rms) be¬ 
tween our results and those of Soubiran et al. (1998) (exclud¬ 
ing HD 175305, for which uncertainties in Soubiran et al. 1998, 
are much larger than typical) and obtained A7’ c n = 63 K, 
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Table 4. Atmospheric parameters taken from Butler et al. (2006) and our results for eight planet-hosting stars. 


Name 

Sp. 

type 

7eff 

[K] 

Butler et al. 

log g 

[cm s~ 2 ] 

[Fe/H] 

T t « 

[K] 

Our results 
iogg 
[cm s -2 ] 

[Fe/H] 

HD 10697 

G5IV 

5680 ± 44 

4.12 ± 0.06 

0.19 ± 0.03 

5611 ± 8 

4.03 ± 0.02 

0.06 ± 0.05 

HD 38529 

G4V 

5697 ± 44 

4.05 ± 0.06 

0.45 ± 0.03 

5620 ± 13 

4.03 ± 0.03 

0.31 ±0.06 

HD 72659 

GO 

5920 ± 44 

4.24 ± 0.06 

-0.01 ± 0.03 

5931 ± 15 

4.28 ± 0.04 

-0.04 ± 0.06 

HD 74156 

GO 

6068 ± 44 

4.26 ± 0.06 

0.13 ±0.03 

6060 ± 23 

4.32 ± 0.05 

0.04 ± 0.05 

HD 75732 

G8V 

5235 ± 44 

4.45 ± 0.06 

0.32 ± 0.03 

5265 ± 15 

4.49 ± 0.05 

0.29 ± 0.07 

HD 88133 

G5 

5494 ± 23 

4.23 ± 0.05 

0.34 ± 0.03 

5397 ± 10 

3.97 ± 0.03 

0.23 ± 0.06 

HD 118203 

K0 

5600 ± 150 

3.87 ± 0.10 

0.10 ± 0.03 

5784 ± 15 

3.94 ± 0.04 

0.09 ± 0.06 

HD 209458 

F8 

6099 ± 44 

4.38 ± 0.06 

0.01 ± 0.03 

6089 ± 20 

4.44 ± 0.05 

-0.05 ± 0.06 



5000 5300 5600 5900 6200 

T e « [K] 



3.8 4.0 4.2 4.4 4.6 

log g 



- 0.1 0.0 0.1 0.2 0.3 0.4 0.5 

[Fe/H] 


Fig. 5. Comparison of the atmospheric parameters for eight known planet-hosting stars. Both relations (top) and differences (bottom) between our 
results and those of Butler et al. (2006) of T e g (a), log g (b), and [Fe/H] (c) are presented. The one-to-one relations are presented as solid lines. 


Alog g = 0.15, A [Fe/H] = 0.02, cr D r eff = 54 K, cr D log g = 0.27, 
and cr D [Fe/H] = 0.25. 

Since we used exactly the same spectra as Soubiran et al. 
(1998), the only sources of the differences in the results and es¬ 
timated uncertainties are the numerical procedures, the selection 
of lines, and the EW measurements. We note that our intrinsic 
uncertainties in T t g, log g are 2.6-3.8 times smaller than the re¬ 
sulting scatter in the results. We also note that our results for T e fl¬ 
are systematically 63 K (3 cr) larger, while our results for log g 
are systematically 0.15 dex (2.1<x) estimated larger. However, 
the [Fe/H] results do not show any systematical effects as com¬ 
pared with Soubiran et al. (1998). The scatter in all results, which 
is comparable to uncertainties as determined from Katz et al. 
(1998) and, especially in [Fe/H], is larger than the average un¬ 
certainty as estimated by Katz et al. (1998), suggests that our 
results are more precise than those of Soubiran et al. (1998). 

5.2. Tests with Butler et al. (2006) data 

We performed another test on the quality and reliability of 
our results using our own HET/HRS spectra of eight known 
hosts of planetary systems: HD 10697, HD 38529, HD 72659, 
HD 74156, HD 75732, HD 88133, HD 118203, and HD 209458, 
previously analyzed by other authors. The sample contained 
mainly dwarfs and sub-giants, with spectral types ranging 
from F8 to K0. Our results were compared to the determina¬ 
tions collected in Butler et al. (2006), who included the nearby 


stars and exoplanets data from the Lick, the Keck, and the 
Anglo-Australian Observatory planet searches. The results of 
our analysis for the eight planet-hosting stars are summarized 
in Table 4, where the atmospheric parameters from Butler et al. 
(2006) are also presented (see Fig. 5). The mean intrinsic uncer¬ 
tainties for the atmospheric parameters we obtained and present 
in Table 4 are crT e ff = 15 K, crlog g - 0.04, and <x[Fe/H] = 0.06, 
while the uncertainties presented by Butler et al. (2006) are 
fTefr = 55 K, crlog g = 0.06, and cr[Fe/H] = 0.03. The mean in¬ 
trinsic uncertainties in T e g as given by the TGVIT are 3.6 times 
smaller than estimated uncertainties given in Butler et al. (2006), 
while our estimated uncertainties in the log g and the [Fe/H] 
are comparable. 

After rejecting HD 118203, for which the uncertainty in T e g 
as given in Butler et al. (2006) is three times larger than aver¬ 
age, we calculated the mean differences and scatter (rms) be¬ 
tween our results and those collected in Butler et al. (2006). We 
obtained A T eS = -31 K, Alog g = -0.02, A [Fe/H] = -0.08, 
croTeff = 49 K, cr D log g - 0.12, and cr D [Fe/H] = 0.05. 

Our intrinsic uncertainties in T e s and log g are 3.3 and three 
times smaller, respectively, than the resulting scatter in the re¬ 
sults. The scatter in [Fe/H] is similar to our uncertainty estimate. 

We note that our results for T e g and [Fe/H] are systemati¬ 
cally smaller by 31 K (2cr) and 0.08 dex (1.3cr), respectively, 
while the results for log g do not show any measurable system¬ 
atical effects as compared with Butler et al. (2006). The scat¬ 
ter among all results is comparable to uncertainties determined 
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Fig. 6. Distributions of atmospheric parameters for the 332 stars with complete spectroscopic analysis: 7’eff, log g, v t , and [Fe/H] (panels a)-d)) 
and their intrinsic uncertainties (panels e)-h)). 


by Butler et al. (2006) and by us, especially in [Fe/H], where 
it is smaller than our estimate. This suggests that the results 
presented in Butler et al. (2006) are of similar quality or more 
accurate than ours. 

In concluding our tests, we note that in general the T e §■, 
log g, and [Fe/H] obtained with the TGVIT and those selected 
from Soubiran et al. (1998) or Butler et al. (2006) agree with 
each other. The uncertainties in the metallicities, estimated in the 
TGVIT from the actual Fe abundance distribution as the standard 
deviation of the mean, are realistic. Our [Fe/H] determinations 
agree with those of Soubiran et al. (1998), but are systematically 
1.3cr lower than those presented in Butler et al. (2006). 

The intrinsic TGVIT uncertainties of T e ff are most certainly 
underestimated. The comparison with results from Soubiran 
et al. (1998) suggests that we should increase our <t7’ c h esti¬ 
mates by a factor of 2.6. Similar analysis of scatter between 
our data and those of Butler et al. (2006) suggests that our un¬ 
certainties are underestimated by a factor of 3.3. If we apply 
such corrected <xc 7’eff = 3.3crT e g estimate, our results agree 
with those of Soubiran et al. (1998) and Butler et al. (2006) 
within lcr c . Similar analysis of scatter leads to the conclusion 
that our crlog g estimates should be increased by a factor of 
three. With such uncertainty estimates our log g determinations 
also agree with those of both Soubiran et al. (1998) and Butler 
et al. (2006) well within crclog g = 3.0<xlog g. Since the micro¬ 
turbulence velocity uncertainty is estimated in the same manner 
as for 7’efj and log g, we expect that these uncertainties are also 
underestimated by a factor of three. 

6. Atmospheric parameters 

The results of our determinations of the atmospheric parameters 
for 348 GK giant stars are presented in Table 5 (Cols. 2-5). In 
all cases, the values of T e g, log g, v t , and [Fe/H] are presented 
together with their intrinsic uncertainties. For 332 objects, we 
obtained consistent solutions for each of the stellar parameters 
in (typically) ten iterations of the procedure described in Sect. 4. 

The spectroscopic analysis revealed that five of our pre¬ 
sumable giants are actually dwarfs (TYC 0435-01209-1, 


TYC 0683-01063-1, TYC 1496-01002-1, TYC 3012-00273- 
1, and TYC 4444-00200-1) with log g > 4.0. On the other 
hand, four objects (TYC 2818-00449-1, TYC 2818-00504-1, 
TYC 2823-01028-1, and TYC 2823-01398-1) recognized in the 
SIMBAD as dwarfs are actually K-type giants. 

The values of 7’eff, log g, l\, and [Fe/H] for 332 stars were 
found to stay generally within the range of the TGVIT model 
grids. However, analysis of 24 stars resulted in effective temper¬ 
atures below 4250 K with other parameters within the TGVIT 
range of applicability. As the Kurucz (1993a) models included 
in the TGVIT reach down to 4000 K, we consider these deter¬ 
minations to be realistic but associated with larger uncertainties, 
as will be illustrated in Sect. 6.1. For 16 stars, the resulting atmo¬ 
spheric parameters were either incoherent or outside the range 
of applicability of the TGVIT. A simplified analysis was per¬ 
formed for these stars, as described in Sect. 6.2, and they will be 
discussed separately. 

The obtained 7’eff ranges between 4055 K and 6239 K 
with a median value at 4736 K. The distribution of 7’eff is 
presented in Fig. 6a. Virtually all our stars have 7’eff be¬ 
tween 4000 K and 5200 K (with the vast majority, ~200, be¬ 
tween 4600 K and 5000 K), making them mostly G8-K2 stars. 
The five dwarfs are among the hottest stars in our sample, 
with 7’eff between 4906 K and 6239 K. The intrinsic uncertainty 
distribution in 7’eff is presented in Fig. 6e. 

The derived log g for our sample stars ranges between 1.39 
and 4.78, with the median of 2.66. The distribution of log g is 
presented in Fig. 6b. The majority of our stars (251) have log g 
of 2.0-3.0, making them generally giants. Small fractions of 
19 bright giants with log g < 2.0 and 3 subgiants with 3.7 < 
log g < 4.0 are present in our sample as well. The intrinsic un¬ 
certainty distribution in log g is presented in Fig. 6f. 

The resulting microturbulence velocity reaches values 
from 0.57 km s _1 to 2.49 km s _1 and has a median value 
of 1.4 km s -1 . Most of our stars have i> t between 1.2 km s _1 
and 1.6 km s 1 (Fig. 6c). The value of u t for five dwarfs in our 
sample is scattered over the range 0.57-1.68 km s -1 . The intrin¬ 
sic uncertainty distribution of u t is presented in Fig. 6g. 
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Fig. 7. Relations between r e u, log g, u t , and [Fe/H] for the 332 stars with complete spectroscopic analysis. The uncertainties in all cases 
are presented. 


The metallicity of stars in our sample, [Fe/H], stays 
within —1.0 to +0.45 range, with a median value of -0.15. 
Figure 6d presents the distribution and shows that our stars are 
generally less metal-abundant than the Sun, with most of them 
having [Fe/H] in the range of -0.5 to 0.0. The intrinsic uncer¬ 
tainty distribution in [Fe/H] is presented in Fig. 6h. 

Figure 7 presents a few relations between purely spectro¬ 
scopic atmospheric parameters presented in Table 5. They illus¬ 
trate a general resemblance of our sample to those studied by 
others (cf. Takeda et al. 2008). In Fig. 7a, one can notice a clear 
dependance of log g on T e g. The log g tends to be higher for hot¬ 
ter stars. For stars with T e g > 5000 K, the dispersion in surface 
gravity increases with effective temperature as well. Noticeable 
relations exist between v t and both log g and T e g (Figs. 7b and d). 
The v t decreases with both log g and Teff. In general, [Fe/H] val¬ 
ues seem to be significantly lower for low log g stars (Fig. 7e) 
in our sample, while they reveal uniform distributions in T e g 
(Fig. 7c) and v t (Fig. 7f). 

6.1. Uncertainty estimates 

The mean intrinsic uncertainties of our determinations were 
found to be crT e g = 13 K, clog g = 0.05, cru t = 0.08 km s _1 , 
and cr[Fe/H] = 0.07. As discussed in Sect. 5, all these uncertain¬ 
ties, except [Fe/H], are numerical, resulting from the iterative 
procedures of the TGVIT. The real uncertainties are factor of 
three larger (since our intrinsic uncertainty scaling may change 
in a comparison with a larger reference sample, we present our 
intrinsic uncertainties instead of the scaled ones). In the case 
of [Fe/H], the numerical uncertainties are the standard deviation 
of the mean Fel and Fell abundances, hence realistic. The dis¬ 
tributions of the intrinsic uncertainties are preseneted in Fig. 6 
(panels e-h). There exists no correlation between the intrin¬ 
sic uncertainties and the parameter values obtained for any of 


the atmospheric parameters. The scatter of the intrinsic uncer¬ 
tainties is also uniformly distributed over the respective param¬ 
eters values, but it is worth noting that for the stars with T e g 
below 4250 K, the intrinsic uncertainties are higher, especially 
in log g and [Fe/H]. In Fig. 8, the intrinsic uncertainties of T e g, 
log g, and [Fe/H] as a function of T e g are presented. 

Several stars from the sample demonstrate uncertainties of 
atmospheric parameters at least two times larger than the mean 
values. In the case of [Fe/H], we noted intrinsic uncertainties 
larger than 0.12 for six stars (TYC 0870-00241-1, TYC 1425- 
01506-1, TYC 3012-01520-1, TYC 3304-00408-1, TYC 4006- 
00980-1, and TYC 4428-01582-1). These result from the highest 
scatter of Fe abundances with respect to the excitation potential 
and EWs in the whole sample. 

Five objects (TYC 0096-00005-1, TYC 1425-01506-1, 
TYC 3011-00791-1, TYC 3304-00408-1, and TYC 3930- 
01790-1) have significantly larger T e g intrinsic uncertain¬ 
ties (>30 K). In the case of log g larger than typical intrinsic 
uncertainty (>0.1) was obtained for nine stars (TYC 0096- 
00005-1, TYC 0276-00507-1, TYC 0870-00241-1, TYC 1425- 
01506-1, TYC 3012-01520-1, TYC 3304-00405-1, TYC 3304- 
00408-1, TYC 3930-01790-1, and TYC 4006-00980-1). The 
precision in v t is the worst (i.e., >0.15 km s 1 ) for eight objects 
(TYC 0683-01063-1, TYC 0870-00084-1, TYC 1425-01506-1, 
TYC 1496-01002-1, TYC 3011-00791-1, TYC 3012-00273-1, 
TYC 3930-01790-1, and TYC 4444-00200-1). 

The reason for such increases in uncertainty is not al¬ 
ways clear. For some stars, it may be associated with low 
effective temperatures close to the TGVIT model grids limit 
of T eS = 4250 K (TYC 0096-0005-1, TYC 1425-01506-1, 
TYC 3304-00405-1, TYC 3304-00408-1, TYC 0276-00507-1, 
TYC 3012-01520-1). In several other cases, the relatively 
low signal-to-noise ratio of our spectra might be the reason 
(TYC 3011-00791-1, TYC 4006-00980-1, TYC 0870-00084-1, 
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As a result, only estimates of masses and radii based on spec- 
trophotometric parallaxes derived here are possible in most 
cases. For such estimates, bolometric luminosities are required 
in addition to our log g, T e g, and metallicities. We calculated the 
intrinsic color index (B-V) o and the bolometric corrections BCy 
for our stars from the empirical calculation of Alonso et al. 
(1999). Using photometry and assuming standard interstellar 
reddening with R — 3.1 (Rieke & Lebofsky 1985), we estimated 
luminosities for 57 stars with usable parallaxes. For 275 ob¬ 
jects with unknown or unreliable HlPPARCOS parallaxes, we ini¬ 
tially assumed My from the empirical calculation of Straizys & 
Kuriliene (1981) and estimated luminosities on that basis. 

Propagation of uncertainty was applied to estimate luminos¬ 
ity uncertainty. The adopted luminosities, together with their un¬ 
certainty estimates, are presented in Table 5 (Col. 11). For the 
57 stars with reliable parallaxes, the values obtained here were 
adopted (the last column in Table 1 identifies the stars with the 
HlPPARCOS parallaxes). For the remaining 275 stars, luminosi¬ 
ties were constrained better via fitting to evolutionary tracks, 
as described in the following section. 


Fig. 8. Relations between the intrinsic uncertainties of crT t i r (top 
panel), crlog g (middle panel), and cr[Fe/H] (bottom panel) vs. T tS for 
332 PTPS stars. 

TYC 1496-01002-1, TYC 3012-00273-1, TYC 4444-00200-1, 
TYC 4428-01582-1). 

It is also possible that the larger-than-expected uncertainties 
may be due to an unresolved stellar companion. This is the case 
with TYC 0870-00241-1, a star with excellent quality spectra 
and all atmospheric parameters well within the TGVIT range: af¬ 
ter several epochs of PTPS RV monitoring, it appeared to have a 
stellar-mass companion (Niedzielski et al., in prep.). 

6.2. Stars with incomplete data 

Sixteen stars of our sample have atmospheric parameters that 
are difficult to derive using the Takeda et al. (2005a) method. 
For these stars, the solutions found with the TGVIT were either 
inconsistent or unrealistic. To include at least some of these stars 
in our analysis, we adopted T e g , log g, and initial luminosity es¬ 
timates from Adamow et al. (in prep.) for all of them. Effective 
temperatures were calculated from the empirical calibration of 
Ramirez & Melendez (2005), based on Tycho and 2MASS pho¬ 
tometry. Rough estimates of log g were obtained using the re¬ 
sults of Bilir et al. (2006) and Gelino et al. (2005). For all these 
stars, the average value of [Fe/H] obtained for the 332 stars with 
complete spectroscopic analysis, [Fe/H] = -0.15, was assumed 
(bottom part of Table 5). 

This approach revealed that three stars were indeed too cool 
for the approach of Takeda et al. (2005a,b). The remaining 13 did 
not show any different characteristics in comparison with the 
332 stars with complete spectroscopic analysis, and the reasons 
why they were not applicable to the TGVIT remain to be clari¬ 
fied. Due to the much larger uncertainty of adopted atmospheric 
parameters for these stars, they will be either discussed sepa¬ 
rately or not considered in the following sections. 

7. Luminosities, masses, ages, and radii 

For 271 stars from our sample, no HlPPARCOS parallaxes 
are available and for another four (TYC 0017-01084-1, 
TYC 3663-00622-1, TYC 4006-00019-1, and TYC 3304- 
00090-1), the HlPPARCOS parallaxes n are unreliable (cr n > n). 


7.1. Stellar masses and luminosities 

The stellar masses M/M Q , the radii R/R 0 , and the estimates 
of ages were obtained by comparing the positions of stars in 
the [log L/Lq, log g, log 7’eff] space with the theoretical evo¬ 
lutionary tracks of Girardi et al. (2000) and Salasnich et al. 
(2000) for a given metallicity (Fig. 10). We used all available 
tracks corresponding to eight metallicity values: [Y = 0.23, Z = 
0.0004], [Y = 0.23,Z = 0.001], [Y = 0.24,Z = 0.004], 
[Y = 0.25, Z = 0.008], [Y = 0.273, Z = 0.019] (the solar com¬ 
position), [ Y = 0.30,Z = 0.03], [Y — 0.32,Z = 0.04], and 
[Y = 0.39,Z = 0.07]. To derive M/M e for every star, the tracks 
for the nearest metallicity were applied. Then the maximum- 
likelihood function defined as 



where g/ hs , q'™ ,d . and cr, denote observed and modeled param¬ 
eters and their uncertainties respectively, was minimized over 
the model parameter space. In our case, the number of param¬ 
eters n was set to three because log L/L Q , log g, and log 7’ c n 
were used. The area over which the parameter space was ex¬ 
plored was defined as ±10cr (intrinsic) in log g and log 7’eff and 
as +3cr in log L/Lq. The resulting run of the x 2 over the model 
parameters was found to be typically relatively flat. Therefore, 
the final stellar mass was calculated as a mean value over an 
extended area of ;y 2 < 3 X^m- Using that procedure, we could 
determine consistent stellar masses and realistically evaluate 
the uncertainties. 

The resulting masses range from 0.6 M/M Q to 3.4 M/Mq 
(see Fig. 9b for a histogram). The majority of stars in our sample 
have masses below 2 M a . We identified, however, a significant 
number (63 or —19%) of stars that fall in the intermediate-mass 
range of 2 M Q < M < 1 M Q . Table 5 (Col. 12) presents the final 
adopted masses. 

The uncertainty of stellar mass obtained by comparison with 
the theoretical evolutionary tracks selected here (Girardi et al. 
2000; Salasnich et al. 2000) depends on how accurately the 
position of a star in the [log L/Lq, log g, log 7’eff] space is 
determined. The effective temperatures and gravities are rela¬ 
tively well determined, and the largest source of uncertainty is 
the parallax, i.e., luminosity. In the case of stars for which the 
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Fig. 9. Distributions of luminosities, masses, and radii (panels a)-c)) as well as their uncertainties (panels d)-f)) for 332 stars with complete 
spectroscopic analysis. 


parallax n was precise enough, choosing the metallicity model 
may cause additional uncertainty since we used the existing 
evolutionary tracks. 

For an average red giant with known parallax, we found that 
the mass may be estimated within 0.2 M 0 in the best cases, while 
the mean uncertainty for the whole sample is <xM/M 0 — 0.3. 
The uncertainties may become even larger, ~ 1 M Q , for very con¬ 
fused evolutionary tracks regions at lower effective temperatures 
(log T eS < 3.65). 

No reliable parallaxes were available for most of our stars. 
As a result, luminosity log L/Lq was constrained better than 
the initial estimate by adopting the luminosity from the fits to 
the evolutionary tracks, corresponding to the determined log g, 
log T e ff, and the stellar mass. In Table 5 (Col. 11), such luminos¬ 
ity estimates are presented for 275 stars. The log L/Lq ranges 
from -0.68 to 2.86, with a peak at 1.6 (Fig. 9a). The mean un¬ 
certainty in log L/Lq was found to be 0.19 in the best cases, but 
the average value for the whole sample is cr log L/Lq = 0.23. 

The H-R diagram for 332 PTPS stars with the complete spec¬ 
troscopic analysis is presented in Fig. 10. 

7.2. Age estimates 

As a large number of the stars presented have roughly solar 
masses, a small uncertainty in mass leads to significant uncer¬ 
tainties in the resulting age. However, for the intermediate-mass 
stars, the precision in mass is worse but the estimated ages are 
more robust. The ages of the majority of our program stars are 
therefore uncertain. The estimated stellar ages are presented in 
Table 5 (Col. 14). A wide range is given for many stars, reflecting 
the complex nature of RGC region and the difficulties in age de¬ 
termination for single field objects. Nevertheless, we found that 
a typical star from our sample is 3-5 Gyr old (mean uncertainty 
in age is 1.1-1.5 Gyr). 


7.3. Stellar radii 

The stellar radii were estimated from either spectroscopic log g 
and the adopted stellar mass or from spectroscopic T e g and the 
adopted luminosity. In every case, a parameter obtained from the 
spectroscopic analysis and one from the model fitting were used 
according to the equations 

V GA Mf) 

R/R Q (g,M)= y r , (2) 

Kq 

or 

R/R 0 (T eff ,L)= ^ ( 77 )’ (3) 

where G is the gravitational constant and M Q , R e , L e , 
and r e fro are the solar values of mass, radius, luminosity, and 
effective temperature. 

Maximum uncertainty was estimated in both approaches 
by applying the logarithmic derivative method. In Fig. 11, 
a comparison of radii derived using both methods is presented. 
Although the radii differ significantly for several stars, the gen¬ 
eral agreement between both R/R e estimates is good. The mean 
absolute difference between radii determined from Eqs. (2) 
and (3) is AR = 0.7 R 0 . The largest deviations exist for the most 
extended stars (cf. Fig. 1 1). They most probably result from in¬ 
accurate luminosity estimates for those luminous stars without 
parallax and from a departure from LTE conditions in stars with 
the most extended outer regions. Due to better consistency with 
the evolutionary status, the radii obtained from Eq. (3) are given 
in Table 5 (Col. 13), along with their uncertainties. 

The radii of stars in our sample range from 0.6 R 0 to 52.1 R 0 
(Fig. 9c), with the majority between 9-11 R 0 . The mean un¬ 
certainty in the derived stellar radii was found to be 0.8 R/Rq 
(Fig. 9f). 
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Fig. 10. H-R diagram for 332 PTPS stars with the complete spectro¬ 
scopic analysis. The theoretical evolutionary tracks (Girardi et al. 2000) 
are presented for stellar masses of 1 -3 M G and several metallicities. The 
dotted, solid, and dashed lines correspond to [Fe/H] = -0.4,0.0,0.3, 
respectively. 



R/Rsun(9,M) 


Fig. 11. Comparison of radii obtained from two sets of stellar parame¬ 
ters: log g, M/M b and T e fF, log L/L e for 332 stars. The solid line corre¬ 
sponds to one-to-one relation. 


7.4. Stars with incomplete data 

For the 16 stars with incomplete atmospheric data discussed 
in Sect. 6.2 and presented in the bottom part of the Table 5, 
we applied essentially the same procedure to obtain luminosi¬ 
ties, masses, ages, and radii. However, as only rough esti¬ 
mates of log g and assumed [Fe/H] were available, the resulting 
masses, ages, and radii are very uncertain. Therefore, we ex¬ 
cluded these 16 stars from our discussion in Sect. 10. 


8. Radial velocities 

A cross-correlation analysis with an artificial template was ap¬ 
plied to obtain absolute radial velocities. To construct the cross¬ 
correlation functions (CCFs), the normalized stellar spectra were 
correlated with a numerical mask consisting of 1 and 0 value 
points (Nowak & Niedzielski 2008). The nonzero points corre¬ 
sponded to the positions of 300 nonblended, isolated stellar ab¬ 
sorption lines at zero velocity. In our RV measurements, only 




oRVccf [km s ’] 


Fig. 12. Histograms of the RVs obtained from the cross-correlation 
function (left panel) and their uncertainties (right panel) for 347 stars. 


the first 17 orders of the “blue” spectra, free from telluric lines, 
were used. The CCF was computed step by step for each veloc¬ 
ity point without merging the orders. For every order, the algo¬ 
rithm selected only lines suitable for a given wavelength range. 
The CCFs from all orders were finally added to construct the 
final CCF for the whole spectrum. The final radial velocities 
were measured by fitting a Gaussian function to the CCF for the 
whole spectrum. The uncertainties of the RVs were computed 
as rms/ V"p7 of the 17 RVs obtained from fitting a Gaussian to 
the CCF for each order separately. The mean standard uncer¬ 
tainty obtained this way was ctRVccf = 0.041 km s _I . In Fig. 12, 
the distribution of resulting RVs and their uncertainties is pre¬ 
sented. The absolute precision of presented RV is lower than 
the (tRVccf- The HRS is neither thermally nor pressure stabi¬ 
lized. As a result, the wavelength scale of our spectra is affected 
by these factors, and the real uncertainties are likely to be of the 
order of 1 km s -1 . 

Table 5 (Cols. 6-8) presents the RVs transformed to the 
barycenter of the Solar System with the algorithm of Stumpff 
( 1980) for all stars in our sample (except for TYC 3304-00479-1, 
due to the lack of suitable spectra), together with the epochs of 
observation (MJD). 


9. Comparison with literature data 

Because the vast majority of our stars were studied here for 
the first time, no comparison with literature data is possible. 
Nevertheless, several individual objects were included in other 
surveys, so we can compare our results with those of other au¬ 
thors, obtained with various methods. 

The spectroscopic determinations of T e g were compared 
with effective temperatures derived by Adamow et al. (in prep.) 
from the empirical calibration of Ramirez & Melendez (2005) 
and the Tycho and the 2MASS photometry, where all our 
stars were included. The comparison of results for 332 stars 
is shown in Fig. 13. No significant systematic shift is present 
in the data, but the scatter of 224 K on average increases to¬ 
wards higher effective temperatures. The mean difference be¬ 
tween our temperature and those of Adamow et al. (in prep.) is 
AT e ff = 48 K; only for four objects is the difference significantly 
larger (>740 K). These stars are TYC 0683-01063-1,TYC 3319- 
00170-1, TYC 3930-00665-1, and TYC 3993-01850-1. They are 
all significantly reddened by the interstellar extinction, which 
probably affected the photometric T e g. 

Our adopted values of stellar radii for 332 stars were com¬ 
pared with estimates obtained from the empirical calibration of 
Alonso et al. (2000), based on the angular diameters attained 
from the infrared flux method and distances computed from the 
HlPPARCOS parallaxes. In most cases, the results of R/R 0 were 
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Fig. 13. Comparison of T eff of the 332 PTPS stars obtained in this study 
and taken from Adamow et al. (in prep.). Both relation {top) and differ¬ 
ences {bottom) are presented. The uncertainties of both determinations 
are shown. The solid line corresponds to one-to-one relation. 



R/Rsun (T e ff.L) 


Fig. 14. Comparison of the R/R e for the 332 PTPS stars derived here 
and calculated from the empirical calibration of Alonso et al. (2000). 
Both relation {top) and differences {bottom) are presented. The uncer¬ 
tainties in our determinations are shown for every target, whereas for 
Alonso et al. (2000) radii a typical uncertainty for 10 R Q star is denoted 
in the upper left corner. The solid line presents the one-to-one relation. 


fairly close to each other and comparable within the estimated 
uncertainties. It is worth mentioning that uncertainties in radii as 
estimated from the calibration of Alonso et al. (2000) are large, 
and only a typical uncertainty for a 10 R Q star is presented in 
Fig. 14. The mean difference between our radii and calibrated 
radii is A{R/R 0 ) = -2.2 with a scatter of 5.8 R 0 . It is interesting 
to note that our estimates give systematically lower values for 
the stars with smaller radii, and systematically larger ones for 
the stars with larger radii. This is most probably caused by the 
departure from the LTE. 

In Fig. 15, we present a collective comparison of our param¬ 
eters for several stars with the data available in the literature. The 
results obtained by other authors agree, in general, with our de¬ 
terminations of r eff , log g, [Fe/H], and R/R,.,. In spite of a very 
low number of stars in common, i.e., five for T e g, four for log g, 
eight for [Fe/H], and five for R/R Q , the agreement is obvious, 


especially in the case of metallicity. It appears from Fig. 15, 
however, that the scatter between the determinations of various 
authors is large, which makes the comparison more difficult. 

In Fig. 16, the RVs obtained using the CCF technique are 
compared with literature data for all stars, for which such de¬ 
terminations were available. In all but three cases, our radial 
velocities agree with those of other authors. Large discrep¬ 
ancies exist in TYC 3012-02518-1, TYC 3930-01790-1. and 
TYC 3304-00090-1. The first two stars are members of binary 
systems (Abt 1981; Meisel 1968), where the latter is in addition 
chromospherically active (e.g., Strassmeier 1994). The reason 
for absolute RV discrepancy for TYC 3304-00090-1 is discussed 
in more detail in Adamow et al. (2012). 

The mean differences between the respective parameters are 
summarized in Table 6. We found that on average our determi¬ 
nations agree with all presented literature results within A7' L . ( r = 
-6 K, Alog g = 0.07, A[Fe/H] = -0.07, A R/R Q = -0.6, and 
ARV = -0.029 km s _1 (with the exception of the three discrepant 
stars). The mean scatter of these comparisons represented by the 
standard deviation is equal to 179 K, 0.76 dex, 0.10 dex, 3.8 R Q , 
and 0.841 km s -1 for T e g, log g, [Fe/H], R, and RV, respectively. 
No systematic effects were found. Our results agree with those 
of other authors within the estimated uncertainties. 


10. Discussion 

Several relations between the integral and the atmospheric pa¬ 
rameters are depicted in Fig. 17. From the mass-luminosity re¬ 
lation (Fig. 17a), we infer that relatively high-mass stars tend 
to be brighter. In the case of solar-mass objects, however, we 
see a large dispersion in log L/Lq. A clear correlation exists 
between log L/Lq and log R/R 0 (Fig. 17b), reflecting the well- 
known relation L ~ R 2 7/f ff . In Fig. 17c, a strong log R/R 0 de- 
pendance on log g, resulting from Eq. (2), is visible. One can also 
see that the hotter stars in our sample tend to be of higher mass 
(Fig. 17d), except for the dwarfs, which reveal similar masses 
for a wide range of effective temperatures (~ 1300 K). These ef¬ 
fects are largely a result of our sample definition, which is well 
illustrated in the log L/L Q vs. log T cl y relation (Fig. 10). In turn, 
masses of our sample stars show quite uniform distribution ver¬ 
sus log g and [Fe/H] (Figs. 17e and f). 

In Fig. 17g, we see that our stars are quite uniformly dis¬ 
tributed over the metallicity vs. radius plane. Figure 17h is sim¬ 
ilar to the mass-luminosity relation (Fig. 17a), probably be¬ 
cause the stellar radii presented in Fig. 17b result from the 
derived luminosities. 

The age-metallicity relation for 332 PTPS stars is shown in 
Fig. 18 (top panel). As we stated in Sect. 7.2, the stellar ages 
are uncertain and in many cases were ambiguously determined. 
The lower limit of log Age is presented if the age range is 
given in Table 5 (Col. 14). However, the well-known tendency 
resulting from the Galactic evolution is maintained. In Fig. 18 
(bottom panel), the relation between stellar age and mass is 
presented (the Pearson correlation coefficient is R — -0.937). 
A quadratic fit can be determined for stars with masses larger 
than the solar mass: M/M Q = (1.07 +0.07)(log Age) 2 -(21.42 + 
1.30)(log Age) + (107.93 ± 6.04). Based Fig. 18, we find that 
the giants in our sample present a wide range of evolutionary 
stages. Relatively older stars (>3.2 Gyr) have simultaneously 
lower masses (~1 M Q ) and reveal a higher spread in metal¬ 
licity (0.26). The more massive giants (>1.5 M Q ), with lower 
spread in metallicity (0.16), are younger (<3.2 Gyr). 

From Fig. 10, one can see that our sample is composed 
mainly of regular giants evolving along the RGB and giants from 
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Fig. 15. Comparison of various parameters for stars with the literature data. The uncertainties are presented if available. The solid lines denote the 
one-to-one relations. 



RV ccf [km s- 1 ] 

Fig. 16. Radial velocities derived from our cross-correlation analysis 
compared with those available in the literature. Relation (top) and dif¬ 
ferences ( bottom ) between the RV results are shown. The uncertain¬ 
ties of the RVs are presented if available. The solid line presents the 
one-to-one relation. 

the RGC. Only very few stars appeared to be dwarfs. We applied 
the definitions of RGC as proposed in Jimenez et al. (1998) and 
Tautvaisiene & Puzeras (2009), i.e., 4700 K < 7' c ir < 5100 K and 
1.5 < log L/Lq < 1.8 to our data. Due to the large uncertainty 
in our luminosities, we extended the log L/Lq range by 0.2 
to [1.3-2.0]. As a result, we found 126 stars to fulfill the criteria 
for the clump giants. These stars constitute only about 38% of 
our sample. The mean absolute magnitudes of these stars in four 
bands (My - (0.985 + 0.499) mag, Mj = (-0.700 + 0.502) mag. 


Mh = (-1.178 ± 0.502) mag, and Mk = (-1.282 + 0.508) mag) 
confirm that they are the clump giants. In Fig. 19, we can see 
our RGC stars together with the rest of the sample in Mk vs. 
T eff and M K vs. [Fe/H] plots. The scatter (defined as the stan¬ 
dard deviation) in the absolute magnitudes is large due to the 
lack of parallaxes, and it is not possible to search for relations 
with either metallicity or temperature. Given the observational 
uncertainties in luminosity and our sample definition, we cannot 
state that the observed relative frequency of RGC stars reflects 
the true one. 

The stellar masses obtained here by fitting observed stellar 
parameters to selected evolutionary tracks are model dependent. 
Since the mass-loss was ignored, they also should be considered 
as upper limits. 

11. Conclusions 

We presented the atmospheric parameters (7’eir, log g, v t , 
and [Fe/H]), luminosities, masses, radii, ages, and absolute ra¬ 
dial velocities for 348 stars from the RGC sample of the PTPS. 
For the vast majority of them, these are the first determinations. 

For 332 stars, the complete spectroscopic analysis resulted 
in precise 7’ c( |, log g, v t , and metallicities. The intrinsic uncer¬ 
tainties in the parameters derived are crT e ff = 13 K, crlog g = 
0.05, crv t = 0.08 km s -1 , and cr[Fe/H] = 0.07 (the uncertain¬ 
ties in effective temperatures, gravities and microturbulence ve¬ 
locities as discussed in Sect. 5 are underestimated by a factor 
of three). In the case of 16 stars, for which either data was in¬ 
complete or our spectroscopic approach failed, effective temper¬ 
atures and log g were estimated from the photometric data. 

In our analysis, five stars from the present sample appeared 
to be dwarfs and 343 giants, 126 of which are located in the 
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Table 6. Mean differences and standard deviations between stellar parameters presented in this work and derived by other authors. 


Reference 

Number 

AZ'gff ± CT 

Alog g±cr 

A [Fe/H] ±cr 

A R ± <x 

ARV ±cr 

Remarks 


of stars 

[K] 

[cm s~ 2 ] 


[Re] 

[km s~‘] 


Adamow et al. (in prep.) 

332 

48 ± 224 

- 

- 

- 

- 

~or~ 

Alonso et al. (2000) 

332 

- 

- 

- 

-2.2 + 5.8 

- 

0) 

Baines et al. (2009) 

1 

-178 

- 

- 

-0.7 

- 


Baines et al. (2010) spectroscopy 

1 

-179 

-0.34 

-0.10 

-5.2 

- 


Baines et al. (2010) interferometry 

1 

-78 

- 

- 

-3.1 

- 


Balachandran et al. (1988) 

1 

- 

- 

-0.32 

- 

- 


Bartkevicius & Sperauskas (1999) 

1 

- 

- 

- 

- 

0.113 


Blackwell & Lynas-Gray (1998) 

1 

22 

- 

- 

- 

- 


Borde et al. (2002) 

1 

-57 

-0.74 

- 

- 

- 


Brown et al. (1989) 

2 

103 ± 129 

0.26 ± 0.42 

-0.14 + 0.10 

- 

- 

(n 

de Medeiros & Mayor (1999) 

4 

- 

- 

- 

- 

0.178 + 0.169 

on 

Duflot et al. (1995) 

5 

- 

- 

- 

- 

-0.405 ±1.810 

On 

Famaey et al. (2005) 

19 

- 

- 

- 

- 

0.163 + 0.834 

On 

Gontcharov (2006) 

21 

- 

- 

- 

- 

0.422 ± 1.475 

On 

Hekker & Melendez (2007) 

1 

-184 

-0.74 

-0.11 

- 

- 


Luck (1991) 

1 

448 

2.14 

-0.01 

- 

- 


Luck & Fleiter (2007) 

1 

-120 

-0.13 

-0.05 

- 

- 


Me William (1990) 

3 

-10 ± 176 

0.05 ± 0.84 

0.06 ± 0.09 

- 

- 

0) 

Nidever et al. (2002) 

1 

- 

- 

- 

- 

-0.647 


Nordgren et al. (1999) 

1 

- 

- 

- 

0.8 

- 


Pasinetti Fracassini et al. (2001) 

4 

- 

- 

- 

1.4 ±6.1 

- 

0) 

Ramirez & Melendez (2005) 

1 

-42 

- 

-0.03 

- 

- 


Soubiran et al. (2008) 

5 

- 

- 

0.00 + 0.12 

- 

- 

0) 

Taylor (1999) 

4 

408 

- 

-0.02 ± 0.06 

- 

- 

0) 

van Belle & von Braun (2009) 

1 

-266 

- 

- 

4.9 

- 



Notes. (t) The value after the plus-minus sign is standard deviation (if number of stars to compare is >1). (,) The difference and scatter value 
obtained, with the exception of discrepant stars (see Sect. 9.1 and Fig. 16 for details). 
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Fig. 17. Relations between log L/L 0 , M/M 0 , log R/R e and the atmospheric parameters for the 332 stars with complete spectroscopic analysis. 
The uncertainties in all cases are presented. 


red clump. Our RGC sample was found to be composed of 
stars with T e g between 4055 K and 6239 K and a median value 
at 4736 K (with majority, ~200 of them, between 4600 K and 
5000 K), generally G8-K2 stars. The log g determined ranges 
between 1.39 and 4.78 with a median of 2.66 (the majority of 
our stars, 251, have log g of 2.0-3.0) making them generally 


giants. A small fraction of 19 bright giants with log g < 2.0 
and three subgiants with 3.7 < log g < 4.0 is present as well. 
The microturbulence velocity, v t , for stars from this sample 
ranges from 0.57 km s _1 to 2.49 km s -1 and has a median 
at 1.4 km s -1 . The metallicity, [Fe/H], of stars in our RGC sam¬ 
ple stays within -1.0 to +0.45 range. With a median of -0.15, 


A91, page 14 of 16 



































P. Zielinski et al.: The Penn State-Torun Centre for Astronomy Planet Search stars. I. 



log Age 

Fig. 18. Relations between [Fe/H] and log Age (top panel ) as well 
as M/M 0 and log Age (bottom panel) for the 332 stars with complete 
spectroscopic analysis. Only the uncertainties in [Fe/H] and M/M 0 are 
presented. The lower limit of log Age is presented if the age range is 
given in Table 5 (Col. 14). 
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Fig. 19. Relations between Mk and T e g (top panel) as well as Mk 
and [Fe/H] (bottom panel) for 327 giants. The RGC stars are presented 
as solid circles, while the rest of stars are shown as open circles. 


our stars are generally less metal abundant than the Sun, and 
most of them have [Fe/H] in the range of -0.5-0.0. 


For all 348 stars, luminosities, masses, radii and ages were 
estimated using the atmospheric parameters presented, archive 
photometric data and the HlPPARCOS parallaxes, when avail¬ 
able. The log L/Lq ranges from -0.68 to 2.86 with the maxi¬ 
mum peak at 1.6. The resulting masses range from 0.6 M/M e 
to 3.4 M/Mq, with the majority of stars having masses be¬ 
low 2 M 0 . We identified, however, 63 (or ~ 19%) stars that fall 
in the intermediate-mass range 2 M Q < M < 7 M Q . The radii 
determined range from 0.6 R e to 52.1 R Q . Most of our stars have 
radii of about 9-11 R e . We found that stars from our sample 
are typically 3-5 Gyr old and have mean uncertainties in age of 
around 1.1 -1.5 Gyr. Average uncertainties in derived parameters 
are cr log L/Lq = 0.23, crM - 0.3 M Q , crR = 0.8 R s . 

We find the precision in stellar parameters derived accept¬ 
able for the main purpose of our project, i.e., a massive planet 
search. The precision of atmospheric parameters allows for more 
sophisticated analysis of individual stars in future. However, 
since the lack of the HlPPARCOS parallaxes was identified as 
the main source of the uncertainty in luminosities, masses, radii, 
and ages for most of our stars, we expect that the Gaia 1 satellite 
(Lindegren et al. 1994; Perryman et al. 1997; Bailer-Jones 2002) 
will enable to constrain the intrinsic parameters better. 
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